Gravitational Wave Emission from the Single-Degenerate Channel of Type la 
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The thermonuclear explosion of a C/O white dwarf as a Type la supernova (SN la) generates 
a kinetic energy comparable to that released by a massive star during a SN II event. Current 
observations and theoretical models have established that SNe la are asymmetric, and therefore- 
like SNe II-potential sources of gravitational wave (GW) radiation. We perform the first detailed 
calculations of the GW emission for a SN la of any type within the single-degenerate channel. 
The gravitationally-confined detonation (GCD) mechanism predicts a strongly-polarized GW burst 
in the frequency band around 1 Hz. Third-generation spaceborne GW observatories currently in 
planning may be able to detect this predicted signal from SNe la at distances up to 1 Mpc. If 
observable, GWs may offer a direct probe into the first few seconds of the SNe la detonation. 

PACS numbers: 04.30. Db 04.40.Dg 97.60.Bw 97.60.-s 



Introduction. — The discovery of the Phillips relation 
[l[ enabled the use of Type la supernova (SN la) as stan- 
dardizable cosmological candles, and ushered in a new 
era of astronomy leading to the discovery of the acceler- 
ation of the universe [HHl- Understanding the evolution 
of SNe la may improve the systematics errors associated 
with the the calibration of the Phillips relation and allow 
precision studies of dark energy [J]. 

The single-degenerate (SD) channel of SNe la consists 
of a C/O white dwarf (WD) accreting mass from a non- 
degenerate main sequence or red giant companion. As 
the WD approaches the Chandrasekhar limit, carbon 
burning is initiated in its convective core 0, Q. After 
a few hundred years of this "simmering" phase, unsta- 
ble thermonuclear burning is expected to ignite at one 
or more off-centered points, giving rise to a buoyantly- 
rising, subsonic deflagration flame bubble. 

Both the deflagration-to-dctonation transition (DDT) 
7 1 and the gravitationally-confined detonation (GCD) 
8 , models for the SN la explosion mechanism begin 
with this stage of off-centered ignition. The predictions of 
the two mechanisms depart by the evolution of the bubble 
as it approaches breakout. The DDT mechanism posits 
that a transition from deflagration to detonation is made 
prior to bubble breakout. While the DDT mechanism 
yields results consistent with observation, it requires that 
the density at which the transition is made to be set as a 
free parameter. In contrast, in the GCD model, the flame 
bubble breaks through the surface of the WD, launching 
ash into a surface flow; see Figure [1] The initial deflagra- 
tion phase is relatively inefficient, typically burning only 
a modest fraction of the mass of the WD. Therefore, in 
the GCD model, the WD remains gravitationally-bound, 
leading to a ram-pressure driven detonation at the point 
opposite of bubble breakout, which subsequently unbinds 
the WD. Significantly, because both the GCD and the 



DDT initiate detonations near the edge of the WD, both 
are intrinsically and strongly asymmetrical. 

Spectropolarimetry measurements of lines of interme- 
diate mass ions, including Si II and Ca II, yield greater 
polarizations in the outer layers of SNe la [lOj, strongly 
suggesting asymmetry in the later phases of burning. 
Furthermore, Maeda et al. ll| demonstrated that an in- 



trinsic asymmetry in la ejecta, viewed from a random 
direction, can account for their spectral evolution diver- 
sity. Taken together, these observations are broadly con- 
sistent with the asymmetries predicted by the DDT and 
the GCD mechanisms of SNe la. The combined intrin- 
sic asymmetry and large kinetic energy of the explosion 
motivates us to consider the nature of gravitational wave 
(GW) radiation from SNe la in this paper. As is the 
case with SNe II, which are known to be candidate GW 
sources in the 10 2 — 10 3 Hz region of LIGO-like instru- 
ments (l3 - [l3 | , the asymmetry of SNe la naturally leads 
to the production of GWs. 

We may simply estimate the GW signal strength of 
a SN la using the Newtonian-quadrupolc approximation 
to the Einstein field equations: h ~ GQ/(c A D). Here 
Q is the quadrupole mass moment, D is the distance 
to the source, and G and c are the universal gravita- 
tional constant, and the speed of light, respectively. We 
can place a robust upper-bound on the SN la GW am- 
plitude by assuming that the la mechanism is highly as- 
pherical, such that the non-spherical kinetic energy is 
Q/4 ~ ££f n - E kin ~ 10 51 ergs fH- With this assump- 
tion, wc estimate an upper-bound of the dimensionlcss 
strain amplitude as 10 -20 at a distance of 10 kpc. 

We next estimate the characteristic GW frequency 
range of a SNe la detonation. The detonation speed 
Vdct within a WD of radius R is set according to the 
Chapman- Jouguet detonation condition [lij], and yields 
a characteristic frequency /d e t ~ v^ e %/(2R). Esti- 
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mating the radius R ~ 2000 km as that of a cold, 
near-Chandrasekhar mass C/O WD, and the Chapman- 
Jouguet speed as v^ ct ~ 10 4 km/s, we establish an 
upper-bound to the frequency of the expected GW sig- 
nal at < 2.5 Hz. This simple estimate implicitly assumes 
no pre-expansion occurs during the initial deflagration 
phase, and corresponds to an overluminous SNe la event 
which will yield ~ 1.2A/ of 56 Ni [T3|. A lower-bound 
to the detonation frequency can be estimated using the 
pre-expansion out to R ~ 4000 km required to produce 
the observed nucleosynthetic yield of ~ O.7M 56 Ni in 
normal brightness la events [18| : /dot ^ 1 Hz. 

If the detonation is sufficiently asymmetric, this GW 
signal may be detectable by proposed third-generation 
GW instruments. We next develop a more refined esti- 
mate by post-processing three-dimensional (3D) simula- 
tions of the GCD mechanism Q . 

Methodology. — We simulate the evolution of a 3D hy- 
drodynamical GCD SNe la in the SD channel [9]. The 
progenitor WD is evolved from the onset of deflagration 
through detonation. The simulations employ the Eu- 
ler equations of inviscid, non-relativistic hydrodynam- 
ics, coupled to Poisson's equation for self-gravity, and 
an advection-diffusion reaction model of the thickened 
combustion front. The non-relativistic dynamical as- 
sumption is well-justified here, because the highest bulk 
fluid speeds achieved during the simulation are of order 
the Chapman- Jouguet speed of the supersonic detona- 
tion front, about ~ 0.03c, and the gravitational redshift 
factor - GM/{c 2 r) ~ 1 • 10~ 6 at the surface of the WD. 

The simulations utilize FLASH a modular, 
component-based application code framework created 
to simulate compressible, reactive astrophysical flows 
[22| . The framework supports a block-structured 
adaptive mesh refinement (AMR) grid, and includes 
a directionally-split inviscid hydrodynamic piecewise- 
parabolic method (PPM) solver [23| ■ A Helmholtz equa- 
tion of state describes the thermodynamic properties of 
the stellar plasma, including contributions from black- 
body radiation, ions, and electrons of an arbitrary de- 
gree of quantum degeneracy 24|. A multipole solver 



method determines the gravitational potential from Pois- 
son's equation. Three scalar progress variables track 
the flames and detonation waves, following the stages 
of carbon burning, quasi-static equilibrium (QSE) relax- 
ation, and nuclear statistical equilibrium (NSE) relax- 
ation [Iim. 

The simulation begins with a C/O WD, with a single 
flame bubble ignited slightly off-center. Previous work on 
the GCD has demonstrated robust detonations, arising 
independent of the resolution in the detonation region, 
and for a wide variety of initial bubble sizes and offsets 
[litniUzi- The simulation we present begins with a 1.37 
Mq C/O WD, with a single pre-ignited flame bubble of 
radius 16 km, and 40 km offset from the center of the 
WD. We use a (49152km) 3 computational mesh that en- 



sures boundary effects have no influence during detona- 
tion expansion. To ensure convergence in the GW signal 
we have performed runs at maximum finest resolutions 
of 8km, I2km, and I6km. 

We employ the standard post-Newtonian expressions 
for the gravitational quadrupole radiation field from a 
3D source. In particular, we calculate the dimension- 
less GW amplitudes h + and h x , which correspond to the 
two independent polarizations in the transverse-traceless 
(TT) gauge [28| . We utilize an expression for the second- 
time derivative of the reduced quadrupole mass moments 
Iij, in which the time-derivatives are reduced to spatial 
derivatives using the mass continuity equation, and the 
hydrodynamical Euler equations [28| . This method is for- 
mally one order higher accuracy than evaluating the time 
derivative of the quadrupole moment twice, as shown in 
29]. 

Results. — We present in Figure [5] the product of the 
dimensionless GW signal strength h and the distance to 
the source D, for observers with line of sight on the equa- 
tor (9 = 7r/2, <fi = 0). The simulation begins at t = with 
the flame bubble ignition. The first ~ 2.5 s after bubble 
ignition correspond to the deflagration phase; see Figure 
[TJ During this phase, the flame bubble rises and breaks 
out of the surface of the WD at ~0.9 s leaving behind 
a "stem" that continues to burn beyond breakout. The 
deflagration phase continues as the ash moves over the 
surface of the WD, and powers a jet that gives rise to a 
detonation front which is initiated at a point opposite of 
bubble breakout at ~ 2.5 s. 

During detonation at ~2.5-3 s, the supersonic deto- 
nation front unbinds the WD, and releases a kinetic en- 
ergy Skin ~ 1.4 • 10 51 ergs and a nuclear energy E DUC ~ 
1.34 • 10 51 ergs. The GW signal is dominated by the 
h + polarization, while the h x polarization deviates only 
slightly from zero due to the growth of non-axisymmetric 
instabilities. The h+D signal is largely shaped by sim- 
ple geometric considerations; it is nearly symmetric in 
time about 2.7 s, when the signal reaches a maximum as 
the detonation front passes over the center of the WD. 
The h+D signal achieves a maximum amplitude of 11.8 
cm before undergoing two phases of decline. During the 
first phase, from ~2.7-3 s, the detonation front propa- 
gates over the second hemisphere of the WD. During the 
second phase, at t >3 s, the detonation has completely 
compressed the WD, and the ejecta freely expand into 
the homologous phase of the explosion. 

An analysis of the GW spectral energy density over all 
angles, figure [2] (inset), shows a characteristic frequency 
near 2 Hz, which is consistent with our estimate. Inte- 
gration yields the total energy released as GW radiation 
to be 7.5 x 10 40 ergs (4.2 x 10" 14 of M c 2 ). A time 
analysis of the radiated GW power and the cumulative 
radiated energy, figure El establishes that all appreciable 
GW energy is released during detonation, by ~ 3 s. 

We conducted a convergence study, keeping the initial 
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conditions and parameters of the runs fixed, and adding 
higher levels of maximum refinement, corresponding to 
finest resolutions of f 6 km, 12 km, and 8 km (figured]). In 
terms of the peak h + D signal amplitude at the equator, 
relative to our highest-resolution 8 km run, our 12 km 
and 16 km runs are converged to within within 1.1% and 
19.1%, respectively. 

In figure [4] we compare the characteristic GW strain 
noise for our calculated GCD SNe la events with some 
third-generation GRW detectors currently in planning, 
including the Big Bang Observer (BBO) |l9j, Deci- 
Hertz Interferometer Gravitational Wave Observatory 
(DECIGO) [H, or the Einstein Telescope (ET) a and 
b[2l| . Detection of nearby extragalactic events out to 1 
Mpc could be possible by BBO if the event happens to 
be oriented such that line of sight is close to the equator. 
A more conservative estimate, based upon randomly- 
oriented events, suggests that extragalactic events would 
be observable by BBO out to about 400 kpc. 

Discussion. — Our single-bubble GCD models produce 
intermediate mass elements at a velocity coordinate ~ 
11,000 km/s, creating a layered structure of IME and 
Fe peak (NSE) products similar to observation 17], LL8 1 . 
However, single-bubble GCD models generally underpro- 
duce intermediate mass elements and overproduce 56 Ni. 
Consequently, the models are generally too luminous 
in comparison to Branch normal la events. The Ni- 
overproduction issue may be resolved by beginning with 
multiple simultaneous ignition points [30] , or through al- 
ternative subgrid models of turbulent nuclear burning, 
which will effectively enhance the turbulent burning rate 
during the deflagration phase, and lead to pre-expansion 
of the WD progenitor. Because the later detonation 
phase establishes the kinetic energy of the explosion in 
the GCD model, we do not expect that enhanced pre- 
expansion will significantly impact the strengths of the 
predicted GW signals presented here. Pre-expansion will, 
however, slightly shift the characteristic GW frequency 
of the GCD mechanism downward by a factor of ~ 2. 

Our current theoretical models predict signal strengths 
which would be detectable for galactic SNe la, which 
are relatively rare, with event rates of one per few hun- 
dred years. However, like early models of SNe II, current 
la models have yet to include rotation, convection, and 
other intrinsically multi-dimensional effects which can be 
expected to yield greater signal strengths than predicted 
here. Simultaneous detection in the optical of nearby ex- 
tragalactic SNe within ~ few Mpc with concurrent mea- 
surements, including non-detections, from future GW ob- 
servatories may ultimately yield useful constraints upon 
these future models. 
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FIG. 1. Time evolution of the early phases of the development of the GCD from breakout through initiation of detonation. 
The ash temperature field is volume rendered, and the surface of the white dwarf is shown in green. The top left-frame shows 
the flame bubble just prior to breakout at t — 0.9 s. The top right shows the bubble just after breakout, depicting the ash 
surface flow over the gravitationally-bound WD. The lower- left frame shows the development of the jet at the point opposite 
to breakout. The final frame shows the initiation of the jet at ~ 2.5 s. Figure reproduced from [9j by permission of the AAS. 
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FIG. 2. h + and h x amplitudes when viewed on the equator (9 — -|, = 0). Inset: Spectral energy density. 




FIG. 3. Power radiated (solid) and integrated power (dashed) of emitted GWs. 
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FIG. 4. Strain noise curves for current proposed detectors sensitive to GCD sources. The filled region represents possible GCD 
event signals in the range between 0.9 Hz and 3.0 Hz, corresponding to significant pre-expansion and overluminous events 
respectively. The region is capped by the characteristic signal strength at 10 kpc for a randomly-oriented source. We also plot 
the maximum GCD amplitudes at 10 kpc and 1 Mpc, which occur when viewed along the equator [6 — n/2). 



